Mounding in Epitaxial Surface Growth 
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For the past two decades the Villain-Lai-Das Sarma equation has served as the theoretical frame- 
work for conserved surface growth processes, such as molecular-beam epitaxy. However some phe- 
nomena, such as mounding, are yet to be fully understood. In the following, we present a systematic 
analysis of the full, original Villain-Lai-Das Sarma equation showing that mound forming terms, 
which should have been included initially on symmetry grounds, are generated under renormalisa- 
tion. A number of widely studied Langevin equations are recovered as limits or trivial fixed points 
of the full theory. 
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Crystal surface growth, in which the interface is driven 
by the deposition of new material from a process such 
as molecular-beam epitaxy (MBE) , has been extensively 
studied [l]-II[ . Although the behaviour of surface growth 
in the absence of conservation generically belongs to the 
Kardar-Parisi-Zhang (KPZ) [5[ universality class, fre- 
quently the growth condition for MBE allow the imposi- 
tion of conservation laws that prohibit KPZ behaviour. 
Typically experimental observations of MBE on semicon- 
ductor wafers are compared with numerical simulations 
run on lattice models, popularly the Wolf- Villain (WV) 
[6fl and the Das Sarma- Tamborenea (DT) models, and 
continuum Langevin equations postulated from physi- 
cal considerations. Analytically, epitaxial surface growth 
is almost exclusively modeled using the Villain- Lai- Das 
Sarma (VLDS) [H, 0| equation. Although successful, the 
theoretical framework seems somewhat incomplete. No 
clear picture has emerged over the alleged exactness of 
scaling relations and the mechanism of mound forma- 
tion is unaccounted for. In the following we explain why 
mounding is, and ought to be, observed either transiently 
or stably in MBE/ VLDS models; mound formation arises 
naturally from a term in the Langevin equation that has 
been variously missed, overlooked or discarded in the lit- 
erature. A resolution to apparent disagreements and mis- 
understandings on the nature of the scaling relations and 
coupling renormalisation is offered. A complete theoret- 
ical picture for the behaviour and the different regimes 
of epitaxial surface growth as originally envisaged two 
decades ago is presented, with the full theory for con- 
served surface growth via ideal MBE given by: 

<9 t 0(x,t) =^2V 2 0-^ 4 V 4 + A 13 V(V0) 3 

+ A 22 V 2 (V0) 2 + kV (V0V 2 0) + n(x, t), (1) 

where the field </)(x, t) is the surface (height) displace- 
ment at x in the co-moving frame atop a d dimen- 
sional substrate at time t. Growth is subject to the 
white noise r\ with the usual correlator {r)(x, t)rj[x' ', t')) = 




FIG. 1: The vertices of the three non-linearities in Eq. {TJ. 
Thick lines with arrows denote the bare propagator (— icj + 
i^k 2 + j/4(k 2 ) 2 ) -1 . Arcs and dashes drawn in narrow lines 
denote inner products of factors of k. The vertex corre- 
sponding to re has all factors of k on one side, in the form 
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2T 2 S d (x -x')S(t-t'). Ideal MBE || was proposed as 
"atomistic stochastic growth without any bulk defects 
or surface overhangs driven by atomic deposition in a 
chemical-bonding environment where surface relaxation 
can occur only through the breaking of bonds." This con- 
strains the surface dynamics to obey mass conservation; 
the deterministic evolution being cast as the divergence 
of some current, thus ruling out a KPZ term. 

Re-writing two of the couplings in a computationally 
convenient form using A 22 = A 22 — re/2 and K = re: 



(A 22 - |) V 2 (v^) 2 + V (V0V 2 0) 

■ k v(v0v 2 0) - \^ 2 {y<t>) 



A 22 V Z [V4> 



(2) 



Giving rise to the vertices shown in Fig.JTJ The original 
formulation of this Langevin equation [J describing 
epitaxial growth on a d dimensional substrate did not 
contain the n term Eq. ([2]). However it should appear 
in two ways. Either it should be there from the start 
from the same symmetry and conservation arguments 
that produce A 22 . Or if it is absent it is generated any- 
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way from A13 and A22 under renormalisation of the full 
theory Eq. (JTJ) . It is this k non-linearity that provides 
a natural mound forming mechanism. The k notation 
has be used in keeping with the term's appearance in an 
unrelated restricted solid-on-solid model [10( ■ 

The original VLDS formulation was derived by con- 
sidering the most general fourth order equation consis- 
tent with the symmetries of the problem, with further 
refinements coming from other physical insights. Later a 
second way of deriving the equation was sought by con- 
sidering the transition rules on a lattice model of the sur- 
face 



dimensionlcss couplings read: 
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14j , finding ways to turn them into a continuum 
equation and dropping terms considered to be irrelevant 
for the scaling behaviour (seeking to capture the "qualita- 
tive features of the surface morphology" [HI)- Virtually 
all of the analysis of the continuum equations has been 
done using the Dynamic Renormalisation Group (DRG), 
a technique that is an extension of Wilson's approach 
to Renormalisation Group (RG) calculations, in that a 
cutoff is used and momentum shells are integrated out 
to find the differential flow equations of the couplings 
(e.g. [H). Instead of using DRG one can instead turn the 
Langevin equation in question into a field theory, the pro- 
gramme pursued below. This entails forming an action 
constrained by the Langevin equation in order to pro- 
duce a (moment) generating functional. This approach 
was developed independently by Janssen [lU and De Do- 
minicis [16j . 

Prosecuting a rcnormaliscd field theory of Eq. (Q]) to 
one loop various theories of epitaxial growth are recov- 
ered as fixed points or limiting cases. In all dimensions 
Eq. ((T|) displays Edwards- Wilkinson (EW) behaviour for 
i>2 greater than the critical point 1/%, with v% = in the 
absence of non-linearities. For vi < the bare prop- 
agator acquires a pole at the characteristic wavelength 
\J V/±/i>2- In dimensions above the critical dimension, 
d c = 4, Eq. ([T]) is trivially governed by Mullins-Herring 
(MH) behaviour at critical vi = v^. However, in dimen- 
sions d — d c — e < d c , the non-linearities A13, A22 and k 
are all (equally) relevant and produce the non-trivial scal- 
ing behaviour characterised below. In the presence of A13, 
the critical becomes negative and there is no longer a 
generic mechanism guaranteeing V2 = as exists in the 
case u% = 0, i.e. the non-trivial scaling is visible only 
after tuning to a critical point. This is the reason why 
A13 has been dismissed originally. However, as pointed 
out by Haselwandter and Vvedensky [13( , the non-trivial 
behaviour of the full theory might still be visible on an 
intermediate scale, beyond which EW might rule. In par- 
ticular, given that k is generated under renormalisation, 
finite mounding is expected to occur generically. What 
is normally referred to as the VLDS equation, however, 
is Eq. ([1]) with V2, A13, k set to zero. 

The ultraviolet is regularised in a perturbation theory 
in small e = d c — d > 0, where the reparamatcrised, 
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X=f- (3) 



The renormalisation of the couplings is determined by ac- 
counting for all logarithmically divergent diagrams con- 
tributing to the proper vertices I^ 1 ' 1 ), F' 1 ' 2 ' and r' 1 ' 3 ', 
the functional derivatives of the Legendre transform F of 
the cumulant generating function, i.e.: 

(k 1 ,..., W „ +m ;^^4,r 2 ,A 13 ,A22,K) = 



n 



n 



r ([4>], [</>]; ^2, ^4, r 2 ,Ai 3 , a 22 ,k) , 

as derived from the one-particle irreducible, amputated 
diagrams contributing to the corresponding correlation 
function. They give rise to the renormalisation of the 
couplings a in the form a R = Z a a. The infrared, on the 
other hand, is regularised by the mass V2 ^ 0, with v% = 
Z2i / 2(i'~ 2 , renormalisation point = 1 and arbitrary 
inverse length /1. 

The infrared stable fixed point is found as a root of the 
set of beta-functions, f3 a = dlna/dln/i|b a rc, where the 
derivative is to be taken for every coupling a at constant 
bare couplings: 
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where we define f3 a = (3 a /a for convenience. The Wilson 
gamma-functions are correspondingly defined as 7 a = 
dlnZa/dln^lbarc: 



72 = 3fl-, 74 = A ( -x - 3x - 1 ) g, 



-X - 3X - 1 



7s = ( 5 - (2 - e/2)A 
"5 



7a = I 4 — A 



722 = ^g, 7 K 



X - 3X - 1 
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5 + - ) 9, 7i3 = 5.g 

X, 



As the noise does not renormalise at any order, 7r = 0. 
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The scaling behaviour of the proper vertices, 
f< n ' m > (k, W ;^^4,r 2 ,A 13 ,A 22 , K ) 

_ ;-f-(d-4)-f (rf+4)+rf+4- 74 (|(m-n)+X 



p(n,m) / * 



. ^-^-2 Ai3r 2 74 -e+ 71 



(4) 



determines the exponents natural to the field theory: 
f("' m ) (k, W ;^^ 4 ,r 2 ,Ai3,A 22 , K ) 

_ ^--f (d-4+?7)-f (d+4+fj)+d+4+S 



r< 



/ - —■ F 2 A 13 ^22 « i 



Where the hat, •, indicates that the Dirac 8 function from 
momentum conservation by translational invariance has 
been divided out, T^ n ^ = <5(ki + . . . + k„ +m )(5(wi + 
. . . + a; n+m )r(™' m ) eliminating one pair of arguments k, 
oj. Normally, growth exponents (a, z and (3 = a/z) char- 
acterise the approach of stationarity from a flat initial 
configuration and the finite size scaling of the roughness. 
In a field theory this is not particularly germane, so the 
exponents are equivalently (but see 17 |) defined on the 
basis of the two point correlation function: 



G 20 (q, W ) =o|q 

-2 

ill 



-(d+z+2a) 



6|q| 



= -f 20 (q, W ) 



2L 2 



| q |74+4 f ii q 



|q|'% 



74+4 



With suitable metric factor a and b and universal scaling 
function Q{x). In the following, we therefore focus on the 
exponents 
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The simultaneous roots of the beta-functions give the 
fixed points of the theory. The infrared stable one at 



X = —2, A = and g = e/5 
*ives 74 = and 72 = 3e/5 and thus 

r n 1 36 



(7) 
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are the exponents of the full VLDS equation Eq. (QJ at 
the critical point vi =v^. These are the same a and z ex- 
ponents as for MH. With A = this implies both A22 and 
k are zero at the fixed point. The renormalisation of the 
full theory Eq. ffl is driven by A13. In agreement with 
previous results [9( with only A22 we find that A22 and k 
do not renormalise themselves at one loop. The contribu- 
tions from diagrams involving only these two couplings 



neatly cancel at one loop but not at higher orders [18j |. 
While k and A22 do not generate each other, they do mix 
under renormalisation. 

Of the two unstable fixed points, the trivial one, g = 0, 
deserves further attention. As g = implies A13 = 
which causes problems with the definition of A, Eq. it 
is not legitimate to naively read all the gamma-functions 
as zero and extract exact MH behaviour. The scaling 
behaviour of this trivial fixed point is normally referred 
to as the VLDS fixed point, observed either by taking the 
limit A13 — >• or by removing it from the initial Langevin 
equation 0, [TJj], possible as A13 is not generated by A22 
or k. Keeping \ an d replacing A by 



\2 

A 22 



(47T) 



(9) 



(10) 



gives /3. = -e - (3 - e/2)[(5/2) X 2 - 3* - lty and 

c 6 1 £ 6 6 

= — , v — — I , z = 4 , a = — , 

3 2 12 3 3 

where z and a are the traditional one loop VLDS expo- 
nents 0| . If only the A22 non-linearity is present, that is 
X = is taken at the trivial fixed point, then Eq. ([TJ be- 
comes the original VLDS formulation. Interestingly the 
corrections to scaling at two-loops predicted by Jansscn 
fl8| on the basis of A22 alone were found to be too small 
in numerical lattice simulations U 21 1 , this may be due 
to a k correction. Renormalisation does however impose 
bounds on or the fixed point of ip would oblige an un- 
physical 1^4. To one loop, in order to get sensible results 

(5/2)x 2 -3x-l needs to be negative, x G [^^, ^jf^], 
sufficiently lar ge k violates this. However the same expo- 
nents emerge [22| if implemented without consideration 
of this. 

We observe that the k coupling on its own is equiva- 
lent to a model proposed by Escudero [23| to reproduce 
VLDS behaviour. The non-linearity proposed in two di- 
mensions, d X x4> d yy 4> — (Sxj,^) 2 , is exactly the vertex pa- 
rameterised by k in Eq. @. Pursuing the calculation 
with only k is especially straight forward, power counting 
reveals that diagrams for its renormalisation constructed 
solely from the k and the noise vertices are always ultra- 
violet finite; in the absence of any other coupling k is not 
rcnormalised at any order. The only renormalisation is 
of the propagator, with one diagram at one-loop order. 
Hence scaling laws are not corrected to any order, yet 
Janssen's general insight is not wrong jT^|; it is a pecu- 
liarity of having only k that leads to non-renormalisation 
of the coupling, as opposed to a neat cancellation (only) 
at one-loop when A22 is also present. Perhaps other anal- 
yses finding exact scaling laws have inadvertently exam- 
ined this case rather than the VLDS equation. It is now 
quite apparent why this model reproduces VLDS expo- 
nents. It is now also apparent that the model's infrared 
stable fixed point is unphysical (to one loop), and its 
behaviour thus not assessable by perturbation theory 
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As has been observed from its two-dimensional form 
[23j . the k non-linearity favours mound formation. In 
addition to Ehrlich-Schwoebel (ES) barriers expressed 
through i>2 0, 13] it provides a natural mechanism at the 
level of the continuum equations for mound formation. 
In the presence of A13, at the critical point, mounding is 
suppressed on the large scale, Eq. 0, yet visible at and 
below (transient) length scales oc n 5 ^ 2e \ Eq. (j4}. Ap- 
pealing to recent work done on coupling flow in models of 
ideal MBE [l3| provides some theoretical justification for 
transient observation of mounding. At the trivial fixed 
point, on the other hand, VLDS scaling applies generi- 
cally as i/| does not suffer an additive rcnormalisation, 
and mounding may be present on all length scales. 

A connection may also be made to the dynamics used 
in lattice models of ideal MBE. There have been several 
investigations into the link between rules of movement in 
lattice models to the terms in continuum equations that 
represent them, usually concentrating on one-dimension 
[24l | . Hugston and Ketterl 25j showed that going from 
lattice rules that seem intuitive, or even computationally 
convenient, to continuum equations is subtle and fraught 
with unintended consequences. Step edge diffusion [26J, 
appearing in two-dimensions, has been proffered as an 
additional mechanism to ES barriers that leads to un- 
stable mounding, for example in the two-dimensional 



WV model [27(1 - Tne lattice rules for the DT model are 



slightly different and result in EW behaviour instead of 
mounding in two-dimensions. Differences in lattice rules 
may well be the distinction between having n or not in 
the continuum equation for a lattice model. 

Observing that the terms paramcterised by A22 and 
k are equivalent in one-dimension offers an explanation 
as to why k was absent in the original derivation of the 
VLDS equation. It was not generated by renormalisation 
since the analysis immediately focused on A13 = to pre- 
vent vi generation. If ever A13 7^ then A22 would be set 
to zero "without loss of generality", as A13 was deemed 
"more relevant" (28|, further analysis halted by the gen- 
eration of i>i . One analysis LJ, |l4j did not commit this 



omission, but unfortunately missed the diagram generat- 
ing k. A similar commentary applies to the use of master 
equations to generate continuum equations from lattice 
rules. The most telling sign is seen when, in develop- 
ing a continuum equation for the WV model, derivatives 
are rearranged in one-dimension and the result carried 
over to higher dimensions where, however, it is invalid 
(e.g. Eq. (23) from Eq. (21) in [lj|). One might wonder 
whether basic errors in multivariate calculus |29| , such as 
assuming wrongly that V(V0) 3 equals 3(V 2 0)(V(/>) 2 [HI 
have thus far concealed the k term. Fourier transforming 
such terms immediately reveals the problem [29! ] . 

In summary we have shown that the original VLDS 
formulation generates a mounding term, n, overlooked 
in previous studies, partly due to non-linearities being 
deliberately omitted, partly the renormalisation schemes 



employed missed the generation of k and partly from an 
apparent misunderstanding of basic multivariate calcu- 
lus. This k term might effectively capture certain lattice 
rules that give rise to mounding in computer simulations; 
inadvertently studying models with only k present may 
lead to concluding scaling relations are in general exact, 
whereas the presence of k with A22 may account for dif- 
ferences in expected scaling corrections. 
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